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I. INTRODUCTION 



With the advent of the Large Hadron Colhder (LHC) one will have an opportunity to 

study proton-proton, proton-nucleus and nucleus-nucleus collisions at the unprecedentedly 

high energies, ^/s = 14, 9 and 6 GeV per nucleon, respectively While the main physics 

drive of the LHC is the search for Higgs boson, supersymmetry and other physics beyond 

the Standard Model, many ideas of the traditional physics of soft and hard hadron-hadron 

interactions can be tested. In particular, one should be able to address the issue of blackening 

of strong interactions at high energies much better than this can be done at the RHIC and 

Tevatron energies. In this work, the term blackening means the approach of a given partial 

wave its limiting value given by unitarity of the scattering operator. We refer to this regime 

n 

the black body limit (BBL). Specifically, the TOTEM collaboration [2] at the LHC intends 
to study the total, elastic and diffractive dissociation proton-proton cross sections at the 
maximal accelerator energy of ^/s = 14 GeV with the aim to test various models, whose 
predictions depend on the way the BBL is implemented. 

It is commonly believed that phenomena associated with high parton densities are more 
pronounced in nuclei than in free nucleons. In this respect, examining the energy and the 
atomic mass number A dependence of total, elastic and diffractive dissociation cross sections 
in hadron-nucleus scattering, one is expected to see an enhancement of the effects related 
to blackening of the proton-proton interaction. 

In this work, we consider total, elastic and diffractive dissociation proton-nucleus cross 
sections. As a starting point, we use the well-established Glauber- Gribov multiple scat- 
tering formalism 0, 4], which is known to work with a few percent accuracy for total and 
elastic hadron-nucleus cross sections. While the Glauber method is essentially based on non- 
relativistic quantum mechanics, which takes into account only elastic intermediate states, 
its generalization by Gribov within the field-theoretical framework also includes inelastic 
(diffractive) intermediate states. The latter is a manifestation of the increase of the coher- 

n 

ence length associated with the given process with energy [5|. A convenient way to model 
essentia, feature of Mgh-ene.gy had.on =cat«„g . by wo.kig, with eiggn^tate. of the 



scattering operator and by introducing cross section fluctuations [g, [Zl, l8|, M, [lO|, lll |. 

The main goal of the present work is to extend a particular model of cross section fluc- 
tuations summarized in |9i] to the RHIC and LHC energies and to make predictions for the 
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total, elastic and diffractive dissociation proton-heavy nucleus cross sections and discuss the 
approach to the black body regime. 

II. HIGH-ENERGY HADRON-NUCLEUS SCATTERING, GLAUBER FORMAL- 
ISM AND CROSS SECTION FLUCTUATIONS 

In order to define and explain the terms "black body (disc) limit", "unitar- 
ity" , "shadowing" and "diffraction", which we extensively use in this work, it is instructive 
to consider a simple example of high-energy scattering on a completely absorbing spherical 
potential with a radius a in non-relativistic quantum mechanics [jj]. Making usual partial 
wave decomposition, one notices that all partial scattering amplitudes with the angular or- 
bital moments / > /max, where /max = pa and p is the projectile momentum, are zero (no 
scattering). On the other hand, for the partial scattering amplitudes with / < /max, scat- 
tering is maximal in the sense that there is no transmitted wave (there is a shadow formed 
right behind the target sphere) and, hence, the scattered wave equals minus the incoming 
wave, i.e. the partial scattering amplitudes are /; = i/{2p) for / < /max- Using the optical 
theorem, one readily finds the total cross section 

CTtot = 27ra^ , (1) 

which is twice as large as the geometric cross section of the target vra^. One can separately 
calculate the elastic cross section with the result o"ei = 7ra^ and, hence, the difference between 
the total and elastic cross sections, the inelastic cross section, is cTinoi = Tra^. 

These classic results can be understood by noticing that the completely absorbing poten- 
tial of radius a serves as a black body obstacle in the way of the incoming plane wave and 
that one deals with diffraction of the plane wave on a black disc. Then in accordance with 
Babinet's principle of wave optics, the intensity of the scattered or diffracted light (which 
is analogous to ad of our quantum mechanical exercise) is equal to the intensity of light 
scattered in diffraction on the circular opening of size a in an opaque screen, which is pro- 
portional to Tra^. At the same time, the intensity of the absorbed light, which is analogous 
to cTineh is also proportional to vra^, which means that a^i = (Xinei = tto^. The considered 
example shows that the formation of a shadow behind the scattering center leads to diffrac- 
tion. If the scattering potential is a black body, scattering is maximal and the elastic cross 
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section (which is, at the same time, the diffractive cross section) equals half the total cross 
section. The latter is twice as large as the geometric transverse cross section of the target 
black disc. A nice discussion of diffraction in high-energy scattering can be found in ISj. 

In order to show that scattering off the black body is indeed maximal, we recall the 
general condition on the partial scattering amplitudes, which is a consequence of unitarity 
of the scattering operator. 



where accounts for inelastic processes; 9 is the scattering angle. Solving Eq. ([2]) for 
Imfi{Q) and choosing the smaller of the two solutions, we obtain 



From this equation, one sees that the maximal value of Jm//(0) is Imf^'^^{9) = 1/(2 p), 
which is exactly the value of the scattering amplitude in the black body scattering problem. 
One can say that the partial scattering amplitudes for / < /max saturate. While in the consid- 
ered simple example blackening of Imfi leads to the energy-independent total cross section, 
it is not the case in a more realistic situation. For instance, our analysis will demonstrate 
that the total proton-nucleus cross section slowly increases with energy regardless that many 
partial waves reach their constant maximal values. 

In a number of models, which discuss saturation in hard processes, one often assumes that 
the total cross section reaches a fixed maximal value or that partial scattering amplitudes 
reach constant values smaller than the maximal 1/(2 p), see e.g. |14i] . 

The choice of the smaller of the two solutions to Eq. ([2]) is a reflection of the fact that in 
hadron-hadron scattering, the imaginary part of the scattering amplitude is driven by the 
inelastic contribution. 

Turning to hadron-nucleus scattering, we notice that while the target nucleus can be 
better approximated by a completely absorbing black disk than the target proton, it is 
still a poor approximation. A better approach was formulated by Glauber [3j. The target 
nucleus is approximated by a static collection of nucleon scatterers so that the phase of the 
elastic scattering amplitude is a sum of the phases accumulated in each projectile- nucleon 
scattering. This means that if we express the elastic hadron-nucleus scattering amplitude 




(2) 



(3) 
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/^(g) in terms of the profile function r^(6), 



fA{q)='£j d%e^1'TA{h), (4) 
then TA{h) can be expressed in terms of the elementary hadron-nucleon profile functions 

r(fe), 

^{b) = -^ I d'qe-'^-'f{q), (5) 



integrated with the nuclear ground state wave function \l/^(ri, r2, • • • , ^a) 

TA{h) = Jd^n d' f2 . . . r4|^A(ri, F2, • • • , ta)]' (^1 -11(1- 1(6 - sl))^ . (6) 

Equations (j4])-([6]) assume that at high energies the small momentum transfer g is perpen- 
dicular to the direction of the beam, i.e. it is a two-dimensional vector. The corresponding 
conjugated variable is the two-dimensional vector of the impact parameter b. In Eq. ((61), 
the vectors Si are the transverse components of the position of the nucleons r^; f{q) is the 
hadron-nucleon scattering amplitude. For sufficiently heavy nuclei {A > 16) it is permissible 
to neglect the nucleon-nucleon correlations in the ground state nuclear wave function, which 
means that each nucleon moves in the nucleus independently, and to write 

i=A 

|^A(fi,f2, . . . ,f^)|2 = [] p^{fi) , (7) 

i=l 

where the nucleon distribution Pa(^ is normalized to unity. The parameterization of pyi(r) 
is detailed in Sect. lIVI Then the nuclear profile function for a heavy nucleus can be presented 
in the following compact form 

r^(6) = 1 - exp (^-A J d^fpA{r)r{b- s)^ . (8) 

The elementary profile function is readily calculated using the standard parameterization 
for the elementary proton-nucleon scattering amplitude 

j(0 = ^P^tot(g)(l-^r/) ^_^(,)^2/2 ^ 
where cxtot is the energy-dependent total cross section; B{s) is the slope of the amplitude; 



1] = Re f{q)/Im f{q). In our numerical analysis, we use 



5(s) = 10.5 + 0.5 In (s/so) GeV''^ 



(10) 
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where sq = 25 GeV; ri = 7i/2 x 0.0808 = 0.127. 

Evaluating T{b — s) using Eq. ([5]) and substituting the result in Eq. ([8]), we obtain the 
Glauber approximation expression for r^(6) 

r^(6) = l-exp(-A/2ortot(s)(l-2r7)T(6)) , (11) 

where 



= / ^'" ^ 2ms) K<\^- + n ■ '1'' 

In the -B(s) limit, the T(6) function takes a more familiar approximate form, T(6) = 



It is interesting to point out that the profile function rA(&) plays the role of the partial 
scattering amplitude and the impact parameter |6| plays the role of the orbital momentum /. 
As a consequence, the unitarity condition is diagonal in |6| and reads (compare to Eq. ([3])) 

2i?er^(6) = |r^(6)|2 + G^°(6). (13) 

The solution to this equation is 



= , (14) 

where rj^ = ImV ^{h) / ReV . The maximal value of ReV Aih) is unity {rjA vanishes in the 
black disc limit), and, therefore, the Glauber approximation expression for TAih) of Eq. ffTTj) 
trivially complies with the unitarity constraint of Eq. (fT3|) . 

The Glauber formalism offers a convenient scheme for the calculation of various observ- 
ables measured in the hadron-nucleus scattering at high-energies such as the total and elastic 
cross sections 

a':^,{s) = 2jdHReTAih), 

a^f{s)= fdH\TA{b)\\ (15) 



It is important to note that while the nuclear profile function saturates, the scattering cross 
sections in Eq. f[T^ grow with energy at large s. 

The quantum mechanical expressions of the Glauber formalism imply that coherent 
diffraction on nuclei consists of only elastic scattering. This contradicts experiments on 
diffraction dissociation, which showed that the incoming particle can dissociate into states 



with the same quantum numbers leaving the target nucleus in its ground state. Therefore, 
the Glauber formalism should be extended to accommodate this experimental fact. 

A simple picture of diffractive dissociation was suggested by Feinberg and Pomer- 
anchuk jsj and elaborated on by Good and Walker One thinks of the incoming wave as 
a coherent superposition of eigenstates of the scattering operator. Each eigenstate interacts 
with the target with its own cross section. Since in general these cross sections (eigenvalues) 
are different, the final state contains not only the initial particle but also other states, which 
diffracted into existence. It is important to note that the formalism of scattering eigenstates 
is based on the assumption that one can represent scattering as superposition of scattering 
of the components with different interaction strengths. The use of this assumption and the 
completeness of the set of scattering states allows to obtain compact formulas. In pertur- 
bative QCD, this assumption can be justified for t ~ relevant for the scattering off nuclei, 
while it is not valid for sufficiently large t. 

Introducing the probability to interact with a given cross section a, P{cr, s), the expres- 
sions for the total and elastic hadron-nucleus cross sections become (compare to Eqs. (fT5l) ) 



a'^^is) = 2 J daP{a) J d'^h ReV Aih.a) 
a,V(s) = j J daP{a)TA{b,a) 



(16) 



In these equations, the profile function T^ib, a) depends on the eigenvalue a rather than on 
the total cross section cr-^otis), 

TA{b, a) = 1 - exp {-A/2 a{l - tr])T{b)) . (17) 

Therefore, the cross sections in Eq. flTB]) are sensitive not only to the first moment of P(o", s), 
(o")(s) = o"f(^(s), but also to higher moments {a^){s). 

The motivation to introduce cross section fiuctuations is the need for a simple picture of 
diffractive dissociation. The cross section for coherent diffraction dissociation of hadrons on 
a nuclear target is found as the difference between the coherent diffraction and elastic cross 
sections |1C| . 



^DDi.s)= d'bl daPia,s) r^(6,a 



2 



daP{a, s)TA{b, a) 



2^ 



(18) 



Since cr'j^{s) is identically zero if cross section fluctuations are absent, cr^£)(s) is the most 
sensitive observable to cross section fluctuations. 
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At small impact parameters and large a, the nuclear profile function saturates, T^ib, a) ~ 
1, and becomes independent of a. This leads to vanishing cr^£)(s). Therefore, cross section 
fluctuations indicate how close to the black body limit regime one is: The proximity to the 
blackening regime is indicated by the decreasing size of cr^{s). Phenomenologically this 
fact can be taken into account by modeling which becomes narrower as ^/s increases 
and by taking into account the increase of cr^^tis) with energy, see Fig. [TJ 

III. ENERGY DEPENDENCE OF P{a, s) 

The distribution over cross sections P{a, s) has the following properties j^: 



The first equation is probability conservation; the second equation requires that P{cr, s) 
reproduces correctly the total hadron-nucleon cross section; the third equation introduces 
ujcr{s) which measures the broadness of cross section fluctuations around the average value. 
One can also consider higher moments of P(cr, s). 

Equations f|T9|) constitute the minimal set of constraints on P(cr, s) and one can success- 
fully model P(cr, s) using only these constraints and the behavior of P(cr, s) in the limiting 
cases of cr — s> and a — > oo. The constituent quark counting rules suggest that -P(cr) = 0{a) 
as 0" ^ 0. In addition, convergence of integrals for the moments of P(cr, s) (see Eqs. (1191) ) 



We assume a particular parameterization of P{cr, s) [9] and determine free parameters 
of the parameterization using Eqs. (fT9|) with crtot(s) and Uo- as an input at each energy. In 
particular, we use the following form for the proton P{a, s), 



whose parameters at typical energies are summarized in Table HI 

It is worth emphasizing that for large and for the nuclear observables considered in 
this paper, effects of fluctuations are primarily determined by the second moment of P(o", s), 
i.e. by the value of the dispersion Ua- [10]. This allows us to use a simple form of P(cr, s) with 




(19) 



requires that P{a, s) ^ faster than any power of a as cr oo. 




(20) 
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Vs, GeV 




Q{s) 


cjo(s), mb 


24 (nD data, [18]) 


0.29 


2.2 


19 


61 (pD data, [19]) 


0.33 


3.4 


16 


546 (UA4, [20]) 


0.19 


0.94 


48 


546 (CDF, [21]) 


0.16 


0.77 


51 


1,800 (CDF, [21]) 


0.15 


0.72 


63 


9,000 (LHC, pA) 


0.10 


0.52 


88 


14,000 (LHC, pp [22]) 


0.065 


0.39 


97.5 



TABLE I: Parameters P{cr, s) at various typical energies. 

energy-dependent parameters, which still captures the essential features of the distribution 
over cross sections. 

The total proton-proton cross section crfotl-s) is calculated using the Regge theory moti- 
vated fit by Donnachie and Landshoff [15-], 



pp I 

tot I 



21.7s 



0.0808 



+ 56.08 s 



-0.4525 



(21) 



which is in a good agreement with the available data. Recently more elaborate parame- 
terizations of the total proton-proton cross section, which explicitly implement Froissart's 
unitarity bound, were suggested 
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171 ]. An inspection shows that all parameterizations 
predict the values of the total proton-proton cross section, which differ by 5-10% at the 
Fermilab and LHC energies. The nuclear cross sections, which we consider, are virtually 
insensitive to such small differences, primarily due to the approximate saturation of the nu- 
clear profile function Fa(&, o"), see the discussion in the end of Sect. HIl We explicitly checked 



that all nuclear cross sections presented in our work change 
the parameterization of cr^otls) of Donnachie and Landshoff 



w at most 1.5%, when instead of 



15| , we use the parameterization 



of [16|]. For the parameterization of [17|, the change is absolutely negligible. 

The parameter uj„ is a key input of our analysis since it defines the broadness of P{a, s) 
[ujcr oc f2(s)) and, hence, the magnitude of cross section fiuctuations. Information on Ua- can 
be extracted either from the inelastic shadowing correction in proton (neutron)-deuterium 
total cross section or from proton-proton or proton-antiproton single diffraction at t = 0, 
see the details in [Ql. For the lower values of ^/s and the UA4 point at ^/s = 546 GeV, 

n 

we used the results of |9[. In particular, there were used the neutron- deuterium total cross 



section data jl8| (with maximal ~ 24 GeV), the analysis of [19| of the proton-deuterium 
data with maximal y/s ~ 61 GeV), and the proton-antiproton single diffraction data taken 
by the UA4 experiment at the SPS collider at CERN with ^/s = 546 GeV Q. 

In addition to this, we used the CDF (Fermilab) data on proton-antiproton single diffrac- 
tion with ^/s = 546 GeV and ^/s = 1800 GeV 2l|]. An extrapolation to the LHC proton- 
proton energy = 14 TeV, Us = 0.06 — 0.07, is done using K. Goulianos fit and is cited 



m 



22l |. A linear interpolation between the y/s = 1.8 TeV and y/s = 14 TeV gives an es- 



timate for the value of Ug at the proton-nucleus LHC energies, ^ 0.10. Note that the 
uncertainty of the extrapolation of diffraction from the Fermilab to the LHC energies (the 
uncertainty in the value of Ug) constitutes the main uncertainty of our predictions for the 
absolute value of c^^)' but it affects only very weakly our predictions for the A-dependence 
of the diffractive cross section. This uncertainty will be rectified during early runs of the 
LHC by the pp experiments which will measure diffraction in pp scattering at small t. 

It is important to note that judging by the the values of Uo- at y/s = 61 GeV and = 546 
GeV, the function Ua- reaches its (broad) maximum around the present RHIC energy range 
of ^/s = 200 GeV. In our analysis, we assumed that ujcr{\^ = 200 GeV) = 0.3. 

Figure [1] shows the distribution P{a, s) as a function of a at three energies considered 
in Table [11 the solid curve corresponds to ^/s = 9 TeV {pA collisions at the LHC); the 
dashed curve corresponds to ^/s = 1.8 TeV (Tevatron); the dot-dashed curve corresponds to 
^/s = 200 GeV (RHIC). As ^/s increases, the position of the maximum of P{cr, s) increases, 
which naturally corresponds to the increasing crfot(s). Although the dispersion becomes 
progressively smaller as the energy increases, there is no significant change in the width of 
the distribution as measured by the range of values of a, where P{a,s) > 0.5 maxP((T, s). 
Consequently, even at the LHC one should expect significant fluctuations in the number of 
wounded nucleons in pA scattering at central impact parameters 23j . 

While the average total cross section increases with energy according to Eq. ( l2T|) . small 
cross sections can grow with ^/s much faster. For instance, the cross sections corresponding 
to P(cr, s) = 0.002 in Fig. [1] increase with energy as cr oc sO-5-o.75_ 
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FIG. 1: The cross section distribution P{a, s) at different energies: the sohd curve corresponds to 
"v/i = 9 TeV (LHC); the dashed curve corresponds to y/s = 1.8 TeV (Tevatron); the dot-dashed 
curve corresponds to ^/s = 200 GeV (RHIC). 



IV. RESULTS AND DISCUSSION 

Using Eqs. (fT3]) and (fTSi) . we calculate the total, elastic and diffractive dissociation cross 
sections for proton-^°^Pb scattering as a function of ^/s. The result is given in Fig. O 

In our numerical analysis, we used the following parameterization of the nucleon distri- 
bution Pa(^) 

P^^'^ = l + exp(Tr-c)/a) ' ^''^ 
where c = Ra — (tt a)^/(3 -R^i) with Ra = 1.145 A^/'^ fm and a = 0.545 fm; the constant po 
is chosen to provide the normalization of pa(^ to unity. 

One sees from Fig. [2] that cross section fluctuations decrease the total and elastic cross 



11 



^ 4000 
g 3800 
o 3600 - 




2000 



1800 



1600 



1400 



1200 



1000 




s^^^ (GeV) 



10^ 10^ lo" 

s^^^ (GeV) 



slOO 




10 



s^'^ fGeV) 



FIG. 2: The proton-Lead total, elastic and diffractive dissociation cross sections as functions of 
y/s. The solid curves correspond to Glauber formalism with cross section fluctuations; the dashed 
curves neglect the cross section fluctuations. 



sections. The effect is largest in the ^/s = 100 — 200 GeV region. This can be explained 
by the increasing role of nuclear shadowing: an increase of c^o- leads to an increase of the 
inelastic shadowing correction, which decreases the total cross section. 

An examination of Fig. [2] shows that, for ^/s > 546 GeV, the total cross section behaves 
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FIG. 3: The atomic number dependence of the total, elastic and diffractive dissociation cross 
sections. The dashed curves correspond to ^/s = 200 GeV and the solid curves correspond to 
^/s = 9 TeV. The second set of dashed and solid curves, which do not go through the points, 
correspond to the approximate calculation of aoD using Eq. (p6]) . 



as 

jA/ \ n.045 



aVtis) oc s'-''' , (23) 

which is slower than the input crff((s) oc s^-'^^'^s. 

The diffractive dissociation cross section (the lower panel of Fig. [2]) noticeably decreases 
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with increasing energies for y/s > 200 GeV. We would like to stress that the predicted 



diffractive dissociation cross section primarily depends on the input Ua- [10| and depends 
only weakly on the shape of the distribution P{a, s). Therefore, the diffractive dissociation 
cross section is a sensitive tool to study the role of cross section fluctuations. 

We also examined the dependence of the total, elastic and diffractive dissociation cross 
sections at ^/s = 200 GeV (RHIC) and ^/s = 9000 GeV {p A at the LHC) on the atomic 
number A. The results are summarized in Fig. [3l where the dashed curves correspond to 
y/s = 200 GeV and the solid curves correspond to ^/s = 9 TeV. 

The total cross section behaves with an increasing atomic number as 

af„^ocA°™ RHIC, 

aliocA'-'' LHC. (24) 

The dependence on the atomic number of the diffractive dissociation cross section is much 
slower 



cr 



ocv4°-27 LHC. (25) 



DD 



The a^jj oc A^'^'' behavior at the LHC kinematics is slower than the oc A^'^ result 



of 10| at much lower energies: cross section fluctuations play a progressively smaller role as 
one increases the energy. 

It was pointed out in [l^ that the fluctuations near the average give the major contri- 
bution to cTfji)- This point was illustrated by Taylor-expanding the integrand in Eq. (|T8ll 
about a = (a) and keeping only first two non-vanishing terms. The approximate expression 
for a'^jj reads Q 

Note that the effects of t] are small and can be neglected. We would like to emphasize that 
the integral in Eq. (l26l) is a smooth function of 6, which does not contain a subtraction of 
two large factors, as appears from Eq. f|T8|) . Therefore, a^^) is much more sensitive to the 
first moments of -P(cr), i.e. to cftotis) and uja{s)i rather than to the details of the shape of 
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Calculations of a^j^ using Eq. fl26|l are presented in the lower panel of Fig. [3]by the second 
set of dashed and solid curves, which do not go through the points. For the LHC energy, 
the approximation of Eq. fl26|) works rather well. For the RHIC energy, the approximation 
of Eq. (!26l) is good only qualitatively. 



V. ELECTROMAGNETIC CONTRIBUTION 



Coherent p A diffraction, p + A ^ X + A, has an important electromagnetic contribution 
originating from the ultraperipheral p A scattering,^ when the nucleus acts as a source of 



quasi- real photons which interact with the proton [24!]. The smallness of the electromagnetic 
coupling constant is compensated by nuclear coherence, which gives the enhancement factor 
Z"^, where Z is the nuclear charge. Therefore, the electromagnetic background becomes 
important for such heavy nuclei as ^osp^ constitutes a correction for light nuclei down 
to 4°Ca. 

Since the strong amplitude is imaginary and the electromagnetic one is real, the two 
contributions do not interfere. Thus, the cross section of this process is given by convolution 
of the flux of the equivalent photons, with the photon-proton cross section, o"^p(co') 

see e.g. 



24| 

<i=/ -nHa7o?M. (27) 



In this equation. 



2Z^a 



M^)^^^ln(^), (28) 

71 \UJ RJ 



where 7 is the Lorentz factor and R is an effective radius of the nucleus; cjmax ^ 1 1 R] ^m\n 
determines the minimal photon energy required to excite an inelastic final state. Assuming 
that the lightest inelastic final state in the 7]? scattering is A(1232), we obtain cUmin = 0.3 
GeV. 

In our numerical analysis of Eq. ( 1271) . for the Lorentz factor in the proton rest frame, we 
used 7 = 2 7l 7l — 1, where 7l and 7l are the laboratory frame Lorentz factors of the nucleus 
and the proton, respectively. This corresponds to 7 ~ 2.3 x 10^ for RHIC and 7 ~ 4.6 x 10^ 
for the LHC. The nuclear effective radius was estimated as i? = Ra = 1.145 A^/^ fm, see 
Eq. (I22I) . The real photon-proton cross section was parameterized in the two-Reggeon form 
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FIG. 4: The electromagnetic contribution evaluated using Eq. (j27p (dashed curves) and coherent 
difFractive dissociation cross sections (solid curves) as functions of ^/s for Pb and Ca. 



of Donnachie and Landshoff [1 



a7„':(s) = 0.0677 s™ + 0.129 s 



-0.4525 



(29) 



where s = 2uj nip + m^. 

The resulting electromagnetic contributions to the coherent diffractive cross section are 
presented in Fig. H] by dashed curves. They should be compared to the coherent diffractive 
dissociation cross sections presented by the solid curves. The comparison shows that the 
electromagnetic contribution completely dominates coherent pA diffraction on Pb-208 at 
all considered energies. For the lighter nucleus of Ca-40, the role of the electromagnetic 
contribution becomes progressively important with an increasing energy: while (T^^ is about 
25% of cr^*^ at the RHIC energy (y^ = 200 GeV), (tI'^ is three times larger than cr^)*^ in 
the LHC kinematics (^i = 9000 GeV). 

VI. CONCLUSIONS AND DISCUSSION 



We calculated the total, elastic and diffractive dissociation proton-nucleus cross sections 
at high energies using the Glauber-Gribov formalism and taking into account inelastic in- 
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termediate states by means of the notion of cross section fluctuations. We extended the 
model of cross section fluctuations of [9] to the RHIC and LHC energies and applied it 
to the calculation of the cross sections. As a consequence of the decrease of cross section 
fluctuations at the LHC energy, we observed a significant reduction of the diffractive disso- 
ciation cross section in pA coherent diffraction. This calculation can serve as a benchmark 
calculation, whose comparison to the future data can give information on blackening of the 
proton-proton interaction. 

We found that towards the LHC energies, y/s = 9 TeV, crf^ oc gO-O'^^, which is slower than 
the input af^t oc s'^-o^'^'^. Studying the cross sections as a function of the atomic number A, 
we found that afft oc A^-™ and ag^ oc A^-^^ = 200 GeV (RHIC) and that afi oc A^-<^^ 

and oc A^ ^^ at = 9 TeV (LHC). 

Another novel result of the present work is an estimate of the electromagnetic contribution 
to coherent p A diffraction coming from ultraperipheral p A scattering. The electromagnetic 
smallness of the background is compensated by nuclear coherence (the enhancement factor 
is proportional to Z^, where Z is the nuclear charge) and the Lorentz 7 factor. We show 
that when the nuclear momentum in the laboratory frame is large, the ultraperipheral 
e.m. background completely dominates coherent pA diffraction on Pb, see Fig. |H One 
way to reduce the electromagnetic contribution is to use lighter nuclei, such as for example 
Ca-40. 
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